Comparison between charge and spin transport in few layer graphene. 
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Transport measurements on few layer graphene (FLG) are important because they interpolate 
between the properties of single layer graphene (SLG) as a true 2-dimensional material and the 
3-dimensional bulk properties of graphite. In this article we present 4-probe local charge transport 
and non-local spin valve and spin precession measurements on lateral spin field-effect transistors 
(FET) on FLG. We study systematically the charge and spin transport properties depending on 
the number of layers and the electrical back gating of the device. We explain the charge transport 
measurements by taking the screening of scattering potentials into account and use the results 
to understand the spin data. The measured samples are between 3 and 20 layers thick, and we 
include in our analysis our earlier results of the measurements on SLG for comparison. In our room 
temperature spin transport measurements we manage to observe spin signals over distances up to 
10 /im and measure enhanced spin-relaxation times with an increasing number of layers, reaching 
Ta ~ 500 ps as a maximum, about 4 times higher than in SLG. The increase of Ta can result from the 
screening of scattering potentials due to additional intrinsic charge carriers in FLG. We calculate 
the density of states (DOS) of FLG using a zone-folding scheme to determine the charge diffusion 
coefficient Dc from the square resistance Rs- The resulting Dc and the spin-diffusion coefficient 
Ds show similar values and depend only weakly on the number of layers and gate induced charge 
carriers. We discuss the implications of this on the identification of the spin-relaxation mechanism. 

PACS numbers: 85.75.-d, 72.80.Vp, 72.25.Dc 



I. INTRODUCTION 



The electronic properties of exfoliated graphene have 
been studied in great detaiPSl^ while the electron spin 
transport still brings up questions. The experimental 
spin-relaxation length As ~ 2 /im in single layer graphene 
(SLG) at room temperature is already promising'^"^ but 
is still at least one order of magnitude below theoretical 
predictiondJISl. As the spin-relaxation is believed to be 
caused mainly by extrinsic scatter ers i n the substrate and 
on the surface of the graphene flak^^, reducing the effect 
of these scatterers should lead to an improvement of the 
electronic transport^ and an increase in A5. 

One way to avoid scatterers is separating the graphene 
flake from the substrate using suspended graphene, re- 
sulting in an increased charge carrier mobility /ii2H12] 
Another way to enhance the transport properties is to 
screen the scattering potentials using few layer graphene 
(FLG). In a stack of graphene layers electrical potentials 
are screened by the o uter la yers with a screening length 
of about 1 to 5 layer^i^li^, depending on the stacking 
ordeJl^. This reduces the effect of external scatterers, 
resulting in only weakly influenced inner layers. The 
effect of screening of the gate induced charge carriers 
on the electrical transport in FLG has b een o bserved 
in several groups' transport measurement^^'^. While 
the spin transport in FLG was examined earliei'^^'^, 
there are no publications on the influence of screening 
on the spin transport properties. The influence of gate 
induced charge carriers on the spin signal is presented 
in Ref . [19| while the effect of those charges on the spin- 
relaxation length still needs to be investigated. Along 



with studying the possible enhancement in spin trans- 
port using FLG, this investigation will also help with un- 
derstanding the effect of possible multilayer inclusions in 
large scale graphene samples in future spintronic devices. 

This paper is organized as follows. In Sec. |ll] we de- 
scribe the selection and preparation of FLG samples. In 
Sec. |III| the charge transport properties of FLG are pre- 
sented. We measure the dependence of the resistance on 
the number of layers and the gate induced charge carri- 
ers and explain the results considering electrical screening 
and a non-uniform background doping of the flake. At 
the end of the section we calculate the density of states 
(DOS) of FLG using a zone-folding scheme. In Sec. IV we 



discuss the spin transport properties of FLG as a func- 
tion of the gate voltage (Vg), compare the results with 
SLG and show the evolution of the spin transport quan- 
tities as a function of the number of layers. Finally we 
compare the spin and charge transport and discuss the 
dominant spin-relaxation mechanism in our devices be- 
fore the paper ends with the conclusions section. 



II. SAMPLE FABRICATION 

The presented measurements were performed at room 
temperature (RT) on mechanically exfoliated FLG flakes 
from highly oriented pyrolytic graphite (HOPG, from 
Advanced Ceramics, AB stacking) on a Si/Si02 substrate 
with an oxide thickness of 300 nm. We determine the 
thickness of the flakes using an atomic force microscope 
in tapping mode (TAFM). The measured thickness t gives 
the number of layers by rounding down the quotient of 
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t and the spacing between two adjacent graphene lay- 
ers dsL = 0.335 nm (corresponding to the thickness of 
SLG), therefore the number of layers is [t/dsLl ■ Due to 
the imprecise nature of the thickness measurements ob- 
tained with TAFM23 and due to the comparison between 
different thickness measurements on the same sample, 
we estimate an error in the number of layers for the FLG 
samples of about 1 layer. Fig.[T](a) shows a scanning elec- 
tron microscope (SEM) picture of a typical sample. The 
illustrated 20-layer FLG flake is contacted with several 
parallel aligned ferromagnetic cobalt electrodes obtained 
with electron beam lithography, e-beam evaporation of 
Co, and a standard lift-off technique. To avoid the con- 
ductivity mismatch and enhance the spin signal, we cover 
the graphene flake with an 0.8 nm thick insulating oxi- 
dized aluminum layer prior to the Co deposition, reach- 
ing contact resistances above Rc = 2 kil. These contact 
resistances are larger than typical FLG resistances on 
a length scale of the spin-relaxation length Ag, a chiev- 
ing in almost all cases non-invasive contactP'^l. The 
highly doped Si-substrate is contacted by an Au elec- 
trode for the electric gating of the device and controlling 
the amount of induced charge carriers Ug in the system. 
The processing is given in detail in Ref. |H 



III. CHARGE TRANSPORT 
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FIG. 1: (Color online) (a) SEM image of a typical spin valve 
device. A 20-layer graphene flake contacted with several par- 
allel ferromagnetic cobalt electrodes (light gray), (b) Nor- 
malized square resistance vs gate voltage for 1, 5, 9, and 20 
layers. 



The samples are first characterized by measuring the 
Vg dependence of the square resistance Rs of the FLG 
flake using local 4-probe geometry. Fig. [l](b) shows three 
typical measurements on FLG and one measurement on 
SLG. All curves show a maximum resistance (minimum 
conductivity (Tmin) at the respective Vg — Vq, marking 
the state where the Fermi energy Ep coincides with the 
lowest DOS. Vg — Vo is therefore the gate voltage with 
the lowest amount of induced charge carriers, correspond- 
ing to Ug = 0. For our samples we get Ug = a{Vg — Vq) 
with a — 7.2 X 10^° cm~^ V~^, calculated using the Si02 
thickness. 

The values of CTmin are displayed in Fig. [2ja) as a 
function of the number of layers. The conductivity in- 
creases approximately linearly with increasing thickness. 
This can also be seen by the fairly constant value for 
Cmin per layer plotted in Fig. [2](b) and can be explained 
by the linear increase of the DOS with the number of 
layers as presented in Fig. [2]jd). The DOS was calcu- 
lated using a zone-folding scheme as described at the 
end of this section and the displayed points show the 
DOS at Ug — 0, corresponding to the minimum value of 
the DOS. The DOS is plotted with and without taking 
into account energy broadening resulting in slightly dif- 
ferent slopes. The linear increase of the DOS points to a 
weak influence of the graphene layer stacking on the DOS 
per layer. For thicker samples we see a small increase 
of (Tmin/layei'2Sl. The conductance per layer of 20-layer 
graphene of fTmin/layer ^ 3.5 x 2e^ /h increases further 
to cr/layer ~ 8.5 x 2e^//i in bulk graphite^^^"^. This rise 



in (Jmin /layer for thicker samples could be explained by 
a stronger influence of the coupling between the layers 
with increasing thickness, which is consistent with the 
interlayer coupling tight binding parameter rising from 
bilayer graphene (BLG) to graphit^^H 

While (Tinin increases with the number of layers, the 
influence of the gate voltage on the resistance is reduced 
with increasing thickness. This can be seen by the in- 
creased full width half maximum (FWHM) of the peak 
shaped resistance curve (see Fig. [ijb) and Fig. [2](c)) 
and can be explained by the distribution of the induced 
charges over the layers. The red, solid curve in Fig. [3](a) 
shows Rs vs Vg measured on a 14-layer sample. Assum- 
ing an equal division of the charges between the SLG-like 
layers, we get a broadened resistance curve (see Fig.|3][a), 
green, dash-dotted curve). Taking the screening of the 
extrinsic potentials (including the gate voltage) by a few 
layer into account, we get even better agreement be- 
tween the modeled resistance curve and the measured one 
(see Fig. |3][a) , black, dotted curve) . The influence of Vg 
on Rs can be described following an easy resistor model 
(see Fig. 4(a) in Ref. [HI). The FLG flake is modeled as 
parallel resistors (the graphene sheets) contacted via an 
interlayer resistance Rint = Pc dsL/-^ at the source and 
drain, where pc ~ 0.1 ficm is the conductance along the 
c-axis of HOPGP^ and A is the contact area. With our 
contact areas of A ~ 0.5 x 0.5 nui^, we get Rint ^ 1 ^, 
much smaller than typical SLG resistances, Rs ^ 2 kfl. 
No further conductance between the layers is considered. 



3 




DOS broadened 
DOS 

10 ' ' 15 ' 20 



300 



200 



100 



300 1<:^ 



200 i> 



100 ■ 



C/5 

o 

Q 



5 10 15 20 5 

# of layers # of layers 

FIG. 2: (Color online) (a) The minimum conductivity, (b) 
the minimum conductivity per layer and (c) the full width 
half maximum (FWHM) of the peak shaped resistance curves 
as a function of the number of layers. The gray points belong 
to a sample that showed overall unusual behavioi'^. (d) The 
DOS at Tig = as a function of the number of layers. The 
blue triangles show the calculated values using a zone-folding 
scheme, and the open black squares show the values including 
an energy broadening of FWIfM « 60 meV. 
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FIG. 3: (Color online) (a) Rs vs Vg for a 14-layer thick sample 
from a 4-probe measurement (red, solid) and from our model 
excluding (green, dash-dotted) and including the screening of 
the gate voltage (black, dotted), (b) Sketch of the model for 
a FLG stack. The blue solid curve shows the modeled SLG 
resistance _Rg°'*"' considering fi — 2520 cm'^V~^s~'^ and a 
minimum carrier density of no = 0.55 x lO^'^ cm~'^. The gray 
dashed curves show the square resistances for the 14-layers 
that, when combined, result in the black dotted curve for the 
modeled stack resistance in (a). 



Using Thomas-Fermi screenin^S^ the total induced 
charge is screened app roxim ately exponentially and dis- 
tributed over the layerfPE^ shown in Fig. ^h). We as- 
sume the resistance per layer to be equal to a modeled 
SLG resistance = (e/x^/n-ind^ + ^^o^)""^, where e 

is the electron charge and n-md and no are the induced 
and the minimum charge carrier densities per layer. The 
total induced charge carrier density rig is distributed 
over the layers, with an exponential decay of the in- 
duced charge carrier density per layer nl^^ (compare 
Ref.HH eq. (5)). Fig. [3];b) shows for SLG with 

/i = 2520 cm^V-^s-i and no = 0.55 x lO^^ cm'^ (blue, 
solid) and the calculated Rs for the different layers in the 
stack (gray, dashed) as a function of Ug. Considering the 
screening length, we follow Ref. [T3]with A = 3 — 5 layers 
and use the best fit to our data, A = 3 layers. In the 
case of SLG it is nind = ng and for FLG X^^ind — '^s- 
The modeled layers closest to the gate still experience a 
strong resistance change by changing n^, while the resis- 
tance of layers farther away is almost unaffected. The 
resulting resistance of the 14-layer stack is plotted in 
Fig. [3]ja) (black, dotted) together with the measured Rs 
of the 14-layer graphene sample (red, solid). We see good 
agreement between the two curves. In the case of mod- 
eling the resistance excluding the screening we have to 
investigate A — ?> oo and get as a result the green dash- 
dotted curve in Fig. [3]ja), which does not fit as well to 
our measurements as the one that includes screening. 



While the shape of the resistance curve for, e.g., the 
14-layer graphene sample can be easily modeled, some 
samples show a further broadened or an asymmetric re- 
sistance curve as a function of Vg (see, e.g., the curve for 
9-layer graphene in Fig. [ijb)). This can be explained as 
follows: Fig.[4]shows the resistance of a 20-layer graphene 
sample, measured on different parts of the FLG flake. 
The black solid curve was measured on a 9 fim long strip, 
while the other three curves represent the resistance of 
sections of this strip. Between the three sections we see a 
shift of Vb by ~ 50 V. Adding up the resistances, results 
for the full distance in the broadened curve (black, solid) 
with a lower maximum resistance compared to a sample 
with a fixed position of Vq for all sections. The shift of 
Vq is caused by a non-uniform background doping of the 
flake that could be due to a locally different doped sub- 
strate, resist residues on the surface of the flake, or the 
metal contacts. This effect could explain the asymmet- 
ric shape of the resistance curves and the spread in the 
values for the minimum conductivity per layer and for 
the FWHM in Fig. ^h) and (c), respectively. The same 
effect has been observed in SLGpSl. 

In the following section we will discuss the spin trans- 
port properties and compare spin with charge diffusion 
in FLG to discuss the spin-relaxation mechanism. There- 
fore, we need to calculate the charge diffusion coefficient 
Dc based on Rs using the Einstein relation 

cr = ei^Dc, (1) 
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FIG. 4: (Color online) 4-probe measurements of the local 
resistance of the 20-layer graphene sample shown in Fig.[l]^a) 
as a function of Vg. The black solid curve shows the resistance 
on a 9 fim long strip (between contact 3 and 7); the other 
curves represent the resistances on parts of this strip, between 
contacts 3 and 5 (5 /im, red, dashed), contacts 5 and 6 (3 fim, 
blue, dotted), and contacts 6 and 7 (1 ^m, magenta, dash- 
dotted) . The vertical lines show the position of the respective 
maximum resistance Vq. 



where iy{E) is the energy dependent DOS and <j{E) = 
1/Rs- In order to easily calculate the DOS for FLG we 
use the fact that the tight-binding Hamiltonian of a FLG 
graphene system can be, in a good approximation, sepa- 
rated into sets of BLG-like and SLG-like Hamiltonian&-°. 
This approach was already experimentally validated by 
infrared spectroscopj^^^^. 

To determine the number and shape of the BLG- and 
SLG-like bands we apply a zone-folding scheme intro- 
duced by Mak et. aP^l that reduces the 3-dimensional 
(3D) band structure of graphite into a 2-dimcnsional 
(2D) band structure for FLG. This approach uses the fact 
that the confinement in the z-direction (perpendicular to 
the FLG flake) induces standing waves and therefore a 
quantization on the wave- vector k^. The quantization 
can be represented by cutting planes in the 3D Brillouin 
zone (BZ) of graphite that cut through different regions 
of the BZ due to the different symmetry-groups for FLG 
with even or odd number of layers^'^. For an odd number 
of layers, there is always a cutting plane through the H 
point, which introduces a linear dispersion band similar 
to SLG. For an even number of layers, such a cutting 
plane is not present. The other cutting planes (if any) 
do not pass through the borders of the graphite BZ and 
introduce BLG-like bands with different effective masses. 

To define the energy dispersion for graphite, we use a 
simple tight-binding approach consisting of only the hop- 
ping parameter for next-neighbors, 70 — 3.15 eV, and an 
interlayer coupling of 71 = 0.37 eV. The inclusion of 
other interlayer and intralayer coupling parameters have 
a minor effect on our results since they are smeared out 
by a broadening introduced by temperature, impurities, 
and other disorder potentials^. To account for such ef- 
fects, we include a Gaussian broadening in the DOS in 
the same way as shown in Ref. 5. The resulting DOS 



at rig = increases linearly with the number of layers 
and is presented in Fig.[2]^d). In the figure the calculated 
values using only the zone-folding scheme are presented 
along with the values including an energy broadening of 
FWHM « 60 meV. This broadening also takes into ac- 
count the effect of the electron and hole pockets around 
Ep present in graphit^^. 



IV. SPIN TRANSPORT 

Now we examine the spin transport properties of 
FLG. Fig. Ista) shows a typical non-local spin valve 
measuremeirP' on FLG. Sending a current / from elec- 
trode 5 to electrode I (see Fig. [TJ^a)) generates a spin 
accumulation at electrode 5. The spins diffuse on both 
sides of the electrode along the flake and generate a volt- 
age drop Vni between electrodes 6 and 9, deflning the 
non-local resistance Rni = Vni/I- Switching the magne- 
tization of one of the inner electrodes (5 or 6) using an 
in-plane magnetic field results in a sign change of R„i (see 
Fig. [sja)). When the outer contacts (1 or 9) are located 
within the spin-relaxation length, additional switches can 
be observed^. The spin valve measurement in Fig.jsjja) is 
taken on a 7-layer graphene sample with an inner contact 
distance of L = 8 /xm. Including the additional switch 
at small field values, we see a spin signal over a distance 
of L = 10 /im. It is worth noting that this is the longest 
distance over which a spin signal has ever been reported 
for graphene based devices. 

For further analysis of the spin transport we perform 
Hanle spin precession measurements^^. They are per- 
formed in the same geometry as the spin valve measure- 
ments with the magnetic field B pointing now perpendic- 
ular to the sample plane causing the injected, in-plane 
oriented spins to precess. The spin dynamics are de- 
scribed by the Bloch equation for the spin accumulation 
,r v24i 



o/o X /tg = 0. 



(2) 



The first term on the left-hand side describes the spin- 
diffusion represented by the spin-diffusion coefficient Dg, 
and the second term describes the spin-relaxation with 
the spin-relaxation time rg. The third term describes the 
precession with the Larmor frequency ujq — giis/fi B, 
where g — 2 is the effective Lande factor and /is is the 
Bohr magneton. In Fig. |5jb) three Hanle measurements 
on a 5-layer graphene sample are presented. Each curve 
consists of the non-local signal acquired for the parallel 
(P) and antiparallel (AP) orientation of the inner con- 
tacts. The black and the blue dots represent the mea- 
surements for L — 2.8 /im and L = 5.4 /im, respectively, 
at the gate voltage Vq- The red curve is measured on the 
longer distance at Vg = Vq — 60 V, where electron charges 
are induced by the gate. The amplitude for the measure- 
ment with increased L is smaller due to additional spin- 
relaxation, as the spins have to travel a longer distance 
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FIG. 5: (Color online) (a) Non-local spin valve signal of a 
7-layer graphene sample. The sweep directions of the mag- 
netic field are indicated (red and black arrows). The distance 
between the inner electrodes is L = 8 ^m. (b) Hanle preces- 
sion measurements of a 5-layer graphene sample at the gate 
voltage resulting in the minimum conductivity Vg — Vo for 
L = 2.8 j-im (black, largest amplitude) and L — 5.4 fim 
(blue, smallest amplitude) and in the hole doped state at 
Vg — Vo ~ 60 V for L = 5.4 /xm (red, intermediate ampli- 
tude). The precession is measured for the parallel (P) and 
antiparallel (AP) configuration of the inner contacts. The 
curve for L = 2.8 fim shows a switch from the P to the AP 
state at —140 mT. This is due to the fact that at relatively 
high fields the non-avoidable in-plane component of the per- 
pendicular field switches the magnetization of one of the inner 
electrodes. 



resulting in a longer time interval for spin-relaxation. In 
addition to the change in the amplitude, a shift in the B- 
field values for the crossing points of the parallel and the 
antiparallel precession curve is visible for the two curves 
measured at Vg = Vq. The crossing points represent the 
-B-field value where the spins have, on average, precessed 
for 90°, resulting in both configurations in a signal of 
Rni ~ 0. An increased distance L, corresponding to an 
increased travel time for the spins, ther efore decreases 
the i3-field which results in 90° -precessiorl^'Sil 

The measurements show that the spin signal can be 
enhanced by inducing more charge carriers (see enhanced 
spin signal comparing the measurement at = Vq — 60 V 
and Vg = Vo). This was also observed for SLGP. 

The Hanle curves can be fitted with the solutions of the 
Bloch-equation ^ , yielding the spin transport quantities 
Ds and rg. Those solutions are calculated with the injec- 
tor considered to be a spin current source and the detec- 
tor considered to be a non-invasive spin voltage probd^*^. 
To exclude (small) spurious background effects, we sub- 
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FIG. 6: (Color online) Ds, ts and Xs (a)-(c) as a function of 
the induced charge carriers Ug for 14-layer graphene and SLG 
and (d)-(f) as a function of the number of layers at = Vq. 
The gray points belong to a sample that showed overall un- 
usual behavioi^^. In addition to Ds, (a) also shows Dc for 
14-layer graphene and SLG and (d) shows the modeled Dc 
at Vg — Vo assuming Rs = 5 kf^/layer. 



tract the AP from the P curve and fit the result. For 
several FLG samples a set of precession measurements 
was performed for different induced charge-carrier densi- 
ties Hg. The spin transport quantities Ds and ts and the 
spin- relaxation length Ag = y/DsTs are plotted as a func- 
tion of in Fig.[6]ja), (b), and (c), respectively. Here the 
results for SLG from Ref. [5]are compared with a 14-layer 
graphene sample representing the results for FLG. The 
general dependence of the quantities is the same for all 
samples. All curves show a minimum at = (corre- 
sponding to = Vo). The change in the three different 
quantities as a function of Ug is minor compared to SLG. 
This can be explained by the fact that in FLG far more 
intrinsic charge carriers are present due to the changed 
band structure compared to SLG, masking the effect of 
the induced charge carriers Ug (see Sec. HI). 

Fig. [6](d), (e), and (f) show Ds, ts, and As, 



respec- 
tively, at Vg = Vq as a function of the number of layers. 
Besides a drop from SLG to BLG, Ds shows no identi- 
fiable dependence and is approximately constant, while 
Ts increases linearly with the number of layers until this 
trend is reduced for the thickest samples, reaching a max- 
imum of Ts ^ 500 ps. With this result, ts still stays far 
below the spin-relaxation times of 20 ns measured at RT 
with electron spin resonance (ESR) in bulk graphit^^ or 
Ts — 55 ns in graphene at T = 150 K^^. There is, at 
this time, no explanation for that difference. The linear 
increase of the measured T5 as a function of the number 
of layers can be explained by the expected screening of 
scattering potentials due to the linear increase of intrin- 
sic charge carriers as a function of the number of layers 
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(compare Fig. [2]^d)). 

The constant value for Dg for more than one layer 
shows that the change in the band structure and the 
screening does not have a strong influence on the spin- 
diffusion. On the other hand the spin-relaxation length 
increases with the number of layers and is doubled at 
Vg = V{) between 1 and 20 layers. The effect of the in- 
duced charge carriers on the spin transport quantities 
is weak for thicker samples. Therefore, we see only a 
small increase of A5 in FLG for rig 7^ compared to SLG 
(Fig. |6][c)) and reach a maximum value of As ^ 3 /im. 
This is also because outer scattering potentials are al- 
ready screened by the intrinsic charge carriers in the 
thicker samples. 

Our results do not show values for A5 as high as re- 
ported by Goto et alW^ (Ag » 8/im). This is probably 
due to the fact that the reported values were derived 
indirectly from spin valve measurements on short dis- 
tances of L ^ 300 nm ^ 8/j,m, making it difficult to 
conclude the behavior over long distances. We also note 
that spin valve measurements are, in general, less conclu- 
sive for spin transport properties than Hanle precession 
measurement^2ll 

Han et al. discuss in Ref. [S] different behaviors of 
the spin signal depending on the induced charge carri- 
ers for different kinds of contact interfaces. We believe 
that in our samples the interface mainly affects the po- 
larization of the injected current and has only weak in- 
fluence on spin-relaxation and spin-diffusion. While we 
have seen reduced spin signals for low Rc in SLG sam- 
ples befor^, the measurements presented here were per- 
formed on samples with Rc values that lead to an R 
parameter of i? > 0.1 /Ltm. As described in Ref. ^ R/\s 
represents the ratio between the contact resistance and 
the graphene resistance over one spin-relaxation length. 
The values that we found for R show that the contacts 
in our samples are non-invasive. Hence, we can rule out 
effects of spins escaping into the cobalt electrodes, fringe 
fields, and interface spin scattering. This is also sup- 
ported by the fact that we have observed spin transport 
under electrically floating cobalt electrodes without any 
measurable effects on the spin signal in SLG and FLG. 

In addition to Ds, Fig. |6][a) and (d) show the charge 
diffusion coefficient Dq, which is calculated using equa- 
tion ^ requiring the DOS. For FLG we use the DOS ob- 
tained by the zone- folding scheme as discussed in Sec. |HI[ 
assuming a broadening of FWHM « 60 meV^^ and, in 
the case of SLG, the broadened DOS discussed in Ref. El 
In Fig.[6]Ja) we use the measured values for the conductiv- 
ity of the samples as a function of V^, while in Fig.lold) we 
assume a fixed resistance per layer of Rs = 5 kI7/layer 
at Vq (compare Fig. ^h)). For SLG the two diflFusion 
coefficients have very similar values independent of Ug. 
In 14-layer graphene Dc is ~ 20% smaller than Ds (and 
^ 50% smaller than in SLG), while both coefficients show 
only a slight change as a function of Ug {Dc oc Dg, see 
Fig. |6][a)). As a function of the number of layers, Dc 
and Ds are approximately constant after the values drop 
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FIG. 7: (Color online) (a) Linear relationship between As 
and Ds, extracted from Fig. [6|a) and (c). The linear fit is 
performed on the red points taken in the electron conduction 
regime, and the black open triangles were taken in the hole 
regime, (b) Slopes of the linear fits of As vs Ds for different 
samples as a function of the number of layers. The increase 
in the slopes is different for the electron and hole conduction 
regimes. The gray point belongs to a sample that showed 
overall unusual behavioi'^. 



between SLG and BLG by roughly ~ 50%. We still see a 
slight decrease from 2 to 5 layers. This shows that here 
Dc behaves the same way as Ds and is not affected by 
the changing band structure or the screening. Similar 
to the case for SLG, the Coulomb electron-electron in- 
teractions still play a minor role in the scattering. The 
main factor limiting diffusion is still impurity potential 
scattering. 

This can also be seen when we plot A5 as a function 
of the corresponding Ds- We observe a linear depen- 
dence for the electron and the hole conduction regimes. 
Fig. ^a.) shows this plot for the 14-layer graphene sample 
(data extracted from Fig.|6]l"^^. Combined with Dc oi Ds 
and A5 = ^/DsTs, this leads to a linear dependence be- 
tween the spin-relaxation time and the momentum scat- 
tering time, which p oints to the Elliott- Yafet-type spin- 
relaxation mechanisnPl2ll, xhe dominance of this mech- 
anism and therefore the spin-relaxation due to impurity 
scattering in our samples can be explained by impurity- 
induced spin-orbit coupling as described in Ref. 1391 
Fig.[7];b) shows \s/Ds for different FLG sampled De- 
pending on whether the points were taken in the electron 
or the hole regime, we see different scalings of the slopes 
with an increasing number of layers. The lower values 
for Xs/Ds in the hole conduction regime indicate that 
there is a higher spin-flip probability for each scatter- 
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ing event. Overall, Xs/Ds increases with the number 
of layers, showing a reduced chance for spin flip with 
an increasing number of layers. This demonstrates the 
enhanced spin transport by using more than one layer 
graphene. 

V. CONCLUSIONS 

We have successfully produced lateral spin field-effect 
transistors (FET) on 3- to 20-layer graphene samples and 
measured the charge and spin transport properties of 
these devices. The reduced influence of external poten- 
tials, including the applied gate voltage, on the charge 
transport with an increasing number of layers has been 
explained by the distribution of the charges between the 
layers. The shape of the resistance curve can be modeled 
if we include screening effects. Further broadening and 
asymmetry of the peak shaped Rs vs Vg curve have been 
explained by inhomogeneous background doping of the 
flakes. Our conductivity per layer for 20-layer graphene 
stays a factor of 2.4 below the value for bulk HOPG. This 
points to extrinsic scattering events at the bottom or the 
top of the flake limiting the transport in our devices. As 
those will be at least partly screened, impurities in the 
layers will also have limiting effects. A weak tempera- 
ture dependence of the resistance (data not shown) also 
points to the fact that impurities and static scatterers 
are the main limiting factors. 

The spin transport quantities Ds and Tg have been 
studied as a function of the induced charge carriers Ug 
and of the number of layers, increases approximately 
linearly with the number of layers showing the expected 
enhancement of the spin lifetime due to the screening 
of the scattering potentials that has been modeled for 
the charge transport. The diffusion coefficients for spin 
{Ds) and charge (Dc) show a decrease from SLG to BLG 
and then stay approximately constant. This shows that 
the number of layers has only a weak influence on the 
diffusion, pointing to a weak coupling between the layers. 



The spin-relaxation length Ag is mainly enhanced for 
Ug = 0. Therefore, we see improvement primarily for 
the spin transport in the regime around amhi- This is 
due to the intrinsic charge carriers of FLG, which mask 
the effect of induced charge carriers. We would like to 
mention at this point that we see no considerable temper- 
ature effect on A5 in FLG (data not shown). This points 
to a negligible effect of phonons on the spin-relaxation. 
The enhancement of As due to screening effects in FLG 
therefore enables the fabrication of improved spin trans- 
port devices. 

Finally, we calculate Dc, using the DOS of FLG ob- 
tained by the zone-folding scheme, and compare the re- 
sult with Dg. As we observe that Dc oc Dg and 
As oc Ds, it seems that the spin-relaxation in our FLG 
samples is mainly due to the EUiott-Yafet mechanism, 
which is also the case for SLG. As the linear dependence 
in our FLG spin transport measurements is based on an 
increase of both values by only a factor of ~ 1.5, this 
result is not yet conclusive and requires further research. 
The theoretical expected dominance of the D'yakonov- 
Perel' spin scattering mechanism^ ^ is probably only mea- 
surable in cleaner sample^^ with higher diffusion coeffi- 
cients and higher mobilities ^. Therefore, measurements 
on high quality graphene spin valve devices have to be 
performed. 

Note added. Related results focusing on the compar- 
ison between spin-relaxation in SLG and BLG have re- 
cently been postecP^l. 
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